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Numerical investigation of the effects of junction geometry on the performance of the fiux-fiow 
type Josephson oscillator has been done. Using direct simulation of the sine-Gordon equation with 
account of RC load the spectral and power properties for different current density distributions have 
been calculated. With discussion of the fundamental differences in the dynamics of the inline and 
overlap junctions the various behavior depending on the junction length and noise intensity has 
been shown. 
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Flux-flow oscillator (FFO) [J, based on a long linear 
Josephson junction, is presently considered as the most 
promising local oscillator for superconducting spectrom- 
eters [2j. For this purpose the spectral linewidth can 
be reduced to 1-10 MHz in the 450-700 GHz frequency 
range. On the other hand, a noisy non-stationary spec- 
trometer Q can be created on the basis of the FFO as 
well. In this case, presently existing FFO designs allow 
reaching a linewidth of more than 50 MHz at generation 
frequency from 450 to 700 GHz. 

The dynamical and spectral properties of the FFO 
were investigated both experimentally [lj, 0,0-0 and 
theoretically IS | . However, experimental and 



theoretical investigations aimed to improving spectral 
and power characteristics were focused on the overlap 
geometry of Josephson tunnel junctions (JTJ) only. Al- 



though in Refs. 1J], 15j the nonuniform current distri- 



butions close to the real experimentally observed profiles 
were studied, orientation of the bias current to the long 
direction of the junction remained intact. The only anal- 
ysis of inline character and comparison for different ge- 
ometries were devoted to the studies of a single fluxon 
motion [lj|-[2l|], escape from the zero- voltage state [22^ . 
Fiske modes [5j and the current-voltage (IV) dependence 
of the FFO regime 23j , 2J] . The oscillation spectra and 
the linewidth of the FFO in the inline geometry has not 
been studied ether experimentally or theoretically. 

The main goal of the present paper is to study the 
influence of the Josephson junction geometry on the flux- 
flow regime and to compare spectral and power properties 
of the inline and overlap bias feeds. 

It is known, that all basic properties of the FFO includ- 
ing spectral characteristics can adequately be described 
in the frame of the sine-Gordon equation: 

(j) tt + acj>t - <j> xx = P4>xxt + v( x ) - sin(0) + r]f(x, t) (1) 

where indices t and x denote temporal and spatial deriva- 
tives, 4> is the phase order parameter. Space and time 
are normalized to the Josephson penetration length Aj 
and to the inverse plasma frequency io~ , respectively, 
a = ujp/uj c is the damping parameter, ui p — y / 2eI c /hC , 
lj c = 2eI c R^ /h, I c is the critical current, C is the junc- 
tion capacitance, Rjy is the normal state resistance, (3 
is the surface loss parameter, n(x) is the injected current 



density, normalized to the critical current density J c , and 
rjf(x,t) is the fluctuational current density. If the critical 
current density is fixed and the fluctuations are treated 
as white Gaussian noise with zero mean, its correlation 
function is: (r]f(x,t)r]f(x',t')) — 2ajS(x — x')8(t — t'), 
where 7 = It/(Jc^j) is the dimensionless noise inten- 
sity, It — 2ekT/h is the thermal current, e is the elec- 
tron charge, Ti is the Planck constant, k is the Boltzmann 
constant and T is the temperature. 

It is well known, that the mode of the current injec- 
tion into the Josephson junctions can affect their static 
and dynamic properties 25]. For the overlap geometry, 
where the bias current is injected into the junction per- 
pendicular to its long direction, most theoretical models 
have assumed uniform distribution of the bias current 
along the junction length. However, in a long narrow 
JTJ the current density distribution across the width of 
a superconducting electrode film is essentially nonuni- 
form. This distribution has singularities at the junction 
edges 25 1, (26|. In the inline case, where the bias current 



is injected into the junction parallel to its long direction 
(which in one-dimensional model leads to injection at the 
ends of the junction only), the dynamic picture is quite 
different. Various junction geometries provide the follow- 
ing current density distributions n(x) [261 ]: 



Vun(x) = ?7o, 



flov(x) = (voL 1 'tt) / 'y 'x(L - x), 
r) m {x) = t] L[5(x) + 5{x - L)], 



(uniform) 
(mixed) (2) 
(inline) 



where i]q is a constant given by the total component of 
the current in the film. The boundaryconditions, that 
simulate simple RC-loads, see Refs. jl3],[i3-[il, have 
the form: 

0x(O,«) + r L c L <f>xt(0,t) - c L <f>tt(0,t) + 

/3r L c L <f> xtt (0 ) t) + (3(t> xt (0,t)=r, (3) 
(f> x (L, t) + r R c R (f) xt (L, t) + c R cf) tt (L, t) + 

l3r R c R xtt (L, t) + pcj> xt (L,t) = T, (4) 

here T is the normalized magnetic field, and L is the 
dimensionless length of the FFO in units of Aj. The di- 
mensionless capacitances and resistances, cl, R and rL,R, 
are the FFO RC-load placed at the left and at the right 
ends, respectively. 



Due to the fact that in the inline case the bias current 
contributes at the ends of the junction only, the inline 
current density distribution of Eq. ([2]) can be accounted 
for in the boundary conditions instead of the sine-Gordon 
equation ((T|), where the term rj(x) is absent (for brevity 
here we do not write the surface losses and RC-load sys- 
tem as it is done in Eqs. © , (0} ) : <p x (0,t) = T- r) L/2, 
(b x {L,t)=T + m L/2. 

At high magnetic fields a dense chain of vortices un- 
der the influence of a bias current is moving unidirection- 
ally through the junction. This fluxon train continuously 
penetrate from one edge of the junction and radiate at 
the other. In the overlap JTJ permanent influence of the 
Lorentz force is maintained through the uniformly dis- 
tributed bias current, while in the inline geometry flux- 
ons receive an accelerating input due to r\i n near the two 
junction ends only. Away from the junction edge this 
inline current influence decreases and for larger damping 
coefficients this effect is reduced more rapidly. As it was 
done for a single fluxon motion [Ii|,[2f|, we can crudely 
estimate the dependence of the FFO regime established 
for a different junction parameters using the following 
terms: in the case aL<Cl the FFO regime in the inline 
JTJ is similar to the overlap case. If aL>l the current- 
voltage dependence as well as spectral properties become 
different for various bias feeds, Eq. ([2]). 
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FIG. 1: FFO linewidth versus total bias current, computed 
for L — 5 and different rj(x), see Eq. ((2|: curve with rectan- 
gles corresponds to uniform distribution of overlap geometry 
r]un, curve with crosses - mixed overlap r\ ov and curve with 
diamonds - inline rji n . Inset: Radiated power from the end x 
= vs 770, the notations are the same as for the linewidth. 

The computer simulations are performed at the same 
parameters, as in work [15| : a=0.033, /3=0.035, T=3.6, 
cl = cr=100, tr—2, rfl=100, but for lower noise inten- 
sity 7=0.05. Firstly let us consider a rather short JTJ 
with L = 5. The IV curves we obtained were similar to 
the results of Ref. Qjjj]. The FFO linewidth for differ- 
ent bias current profiles is presented in Fig. [T] while the 
power at RC load - at the inset of Fig. [TJ As discussed 
above, for inline junctions aL < 1 the vortices starting 
to move from the input edge under the boundary influ- 



ence of a bias current, do not have a large speed reduc- 
tion in comparison with overlap JTJ, so FFO regimes 
established are the same with similar spectral and power 
characteristics. 
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FIG. 2: The IVC of loaded FFO. Uniform overlap (curves 
coincide): L = 20 (curve with crosses) and L = 40 (curve with 
rectangles), mixed overlap: L — 20 (curve with triangles) and 
L = 40 (curve with diamonds), inline: L = 20 (curve with 
stars) and L = 40 (curve with circles). 

Now let us consider more interesting from a practi- 
cal point of view case of a long junction aL > 1. The 
IV characteristics for different junction geometries and 
L = 20, L = 40 are presented in Fig. [2] As it has been 
shown in Ref. [28], for uniform current distribution the 
dependence of the flux-flow step height 770^ on the nor- 
malized junction length L is nonlinear: for small lengths 
the step height 770/1 increases with L and then becomes 
almost constant. As one can see from Fig. [21 the IV 
curves for current density rj un actually coincide. While 
for mixed overlap and inline cases the step heights for 
L = 20 are larger than for L = 40, i.e. from a certain 
L values increase of junction length leads to decrease of 
step heights. In a very long inline junctions the flux-flow 
step may disappear entirely. This is because the energy 
input from such long-distance edges does not allow to 
create a regime of vortex train flow. 

The maximum height of the current step 77071 corre- 
sponds to voltage Vh which can be determined from the 
velocity- matching condition ^28j, so that Vh = T. For 
different distributions ?7 un , ?7 „ and different lengths com- 
puter simulations agrees well with the obtained analyt- 
ical relation, see Fig. O But, in the inline character 
the flow velocity varies spatially [3], and the longer the 
junction, the greater this change. Therefore, even for 
different lengths of inline junctions the voltages Vh are 
different. 

Now let us look at the spectral and power characteris- 
tics for junctions with L = 40. The largest signal in the 
flux- flow regime is obtained biasing the junction near the 
top of the flux- flow branch rjoh. Although the step height 
for various junction geometries is different, the maximum 
radiated power for all three cases is almost the same and 
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FIG. 3: The FFO linewidth for L = 40 and different junc- 
tion geometries: uniform overlap (violet curve and rectangles 
- simulations and theory (0), mixed overlap (red curve and 
crosses - simulations and theory ©), inline (blue curve and 
diamonds - simulations and theory (0)). The inset: the radi- 
ated power (simulations) for the same geometries. 



corresponds to their own 770/1, see the inset of Fig. [3] 
The linewidths of the FFOs at the largest signal are: 
Af un = 0.0179, Af ov = 0.0133, Af m = 0.0247, while 
the minimal attainable linewidths are Af un = 0.0094, 
Af ov = 0.0133, Af in = 0.0243. So, for the inline JTJ 
the linewidth is 2-2.5 times larger than for the overlap 
uniform and mixed current distributions. Therefore, the 
inline junctions seem to be more suitable for an applica- 
tion of a noisy FFO with a broad linewidth and quasi- 
monochromatic generation. However, special designs of 
the overlap junction electrodes can possibly lead to larger 
linewidths than the inline junctions have [8(, but this sub- 
ject is out of scope of the present paper. 

In Ref. [13], the formula for the FFO linewidth, that 
takes into account not only conventional differential resis- 
tance over bias current rj, = dv/dr/o, but also differential 
resistance over magnetic field r^ L = Ldv/dT (control line 
current) was derived: 



Af FFO = 2 ai {r d + ar^ L Y/L 



(5) 



and demonstrated good agreement with experiment 
0,0, and also with computer simulation results for over- 
lap junctions [III], 16]. In Fig. [3] the results of com- 
puter simulations of the linewidth for different bias feeds 
are compared with the theory ^ with the only fitting 
parameter a (violet curve and rectangles - simulations 
and theory ([5]) with a — 0.205 for uniform overlap bias, 
red curve and crosses - simulations and theory ([5]) with 
(7 = 0.06 for mixed overlap bias, and blue curve and di- 
amonds - simulations and theory ([5]) with a = 0.01 for 
inline bias). One can see that theory ^ demonstrates 
good agreement with the computer simulation results for 
all three cases, and that the appropriate value of a de- 
creases with increase of the linewidth, so a is minimal for 
the inline case. Certain disagreement in the area of small 
and large t\q is due to the fact that away from the flux 



flow step the IV curves for different values of T actually 
coincide, and it is impossible to correctly calculate r% L 
in this range. 




FIG. 4: The IVC of loaded FFO, inline geometry, L = 40. 
Curve with stars - T = 1.5, curve with diamonds - F = 2.2, 
curve with rectangles - F — 3.2, curve with crosses - F = 4.2, 
curve with circles - F = 5.2. 
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FIG. 5: Linewidth and radiated power (the inset) of the in- 
line FFO, L = 40. Curve with stars - T = 1.5, curve with 
diamonds - V = 2.2, curve with rectangles - T = 3.2, curve 
with crosses - V = 4.2, curve with circles - T = 5.2. 

It is well known , , 0] , [3 that in the overlap JTJ 
reducing of the magnetic field (T < 2.5) leads to trans- 
formation of the flux-flow regime into Fiske steps, where 
is no possibility of continuous frequency tuning. In Fig. 
[4] the IV curves for inline junction L — 40 depending on 
the external magnetic field T are presented. The FFO 
regime is stored in the broad range T > 1.5 (the Fiske 
steps are almost invisible even in the limit of vanishing 
noise intensity and continuous frequency tuning seems 
to be possible) and the minimal attainable linewidth for 
each r in the working areas is roughly the same, see Fig. 
[5l The maximum radiated power slightly varies versus T 
and the largest signal corresponds to T rs 3.2. 

The influence of noise on the FFO spectral properties 
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FIG. 6: Minimal attainable linewidth and maximal radiated 
power (the inset) for L = 40, V = 3.6 and different junction 
geometries: uniform overlap (curve with triangles), mixed 
overlap (curve with diamonds), inline (curve with crosses). 
Dashed lines - linear fit to computed data. 

is estimated by plotting the minimal attainable linewidth 
and the maximal power values for different current den- 
sity distributions, see Fig. [T5J It is seen that the ef- 
fect of noise on inline JTJ in comparison with overlap 
character results both in the larger linewidth and in the 
steeper slope of the Afmini'y) curve. For the uniform 
and mixed overlap cases the slopes are O.I87 and 0.257, 
respectively, while for inline junction it corresponds to 
O.57 dependence. Also, one can see the deviation from 
linear dependence of the linewidth starting from 7 > 0.1. 
Behavior of the radiated power is strictly non-linear: for 
all geometries it sharply decreases by a factor of three for 
7 = 0.1 compared with 7 = 0.01. 

In conclusions, the spectral and power properties of a 
long Josephson tunnel junction have been studied for var- 
ious junction geometries by numerical solution of the one- 
dimensional sine-Gordon equation taking into account 
surface losses and RC load. The fundamental differences 
in the dynamics of the FFO with various current den- 
sity distributions have been discussed and the distinct 
behavior of inline and overlap junctions with changing 
the external magnetic field and noise intensity is shown. 
Lower performance of the inline junction (compared to 
the overlap) as a local oscillator is observed. On the other 
hand the broad range of magnetic fields, where contin- 
uous frequency tuning is possible, and a large linewidth 
allows one to consider inline junctions as an application 
for a noisy non-stationary FFO spectrometer. 
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